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Abstract 

We report evidence for the B meson decays, B ± J/ipr,^ and B° J/i>r]K Q s , using 90 mil- 
lion BB events collected at the T(4S) resonance with the BABAR detector at the PEP-II e + e~ 
asymmetric-energy storage ring. We obtain preliminary branching fractions in the charged and 
neutral channels of (10.8 ± 2.3(stat.) ± 2A(syst.)) x 10~ 5 and (8.4 ± 2.6(stat.) ± 2.7(syst.)) x 10~ 5 , 
respectively. 
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The first observation of color-suppressed B decay modes with hidden strangeness, ss, was in 
the decay mode, B — > J /ift(ftK, by the CLEO collaboration and more recently from the BaBar 
collaboration '%. The color suppressed B decay modes are represented at the parton level in 
Figure [IJ 



(a) 




(b) 
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Figure 1: The Feynman diagrams of color-suppressed B decays 



Figures^a) and (b) show Feynman diagrams for two-body final state color-suppressed B decays 
into Cabbibo allowed (J/iftK) and Cabbibo suppressed (J/iftir) modes. The color suppression is due 
to the requirement that the color of the quark and anti-quark from the W decay must appropriately 
be matched with the charm quark and the anti-d quark such that the final state particles are color 
singlets. Producing a (ft meson in a three-body final state may be achieved by the addition of ss 
quark-antiquark pairs as shown in Figure ^ ( c )- It is possible also that a (ft is formed from three 
gluon emission as shown in the OZI violating decay in Figure ^d). Since the rj meson is described 
by the quark-antiquark combination, 



V 



uu + dd — 2ss 
72 



we would also expect B mesons, via the similar diagrams, to decay into states with an 77, such as 
B — > J/vfjrjK, where the ss pairs or two gluons form an rj instead of a (ft. Experimentally, the ry is 
seen in charmless decays at a somewhat anomalously high rate which suggests that gluonic 77 
couplings are large. This may lead to an enhancement in the decay rate of B — > J/vftrjK. Since the 
basic diagrams are the same, we would expect the branching fraction of this mode to be comparable 
to BR(B -> J/tft(ftK)= (4.4 ± 1.4 ± 0.5) x 10" 5 2_. If the rate for non-resonant B J/iftrjK is 
found to be much higher or lower, then some other physics may be indicated. 

If the mode B — > J/iftrjK is observed, the neutral mode may be useful for CP violation stud- 
ies jl], j^j and investigation of quasi-two body decays. The decay, B — > ift(2s)K, has been observed, 
and hence the decay B — > ift(2s)K — > J/vftijK should be produced where the J/ipr] are in a rel- 
ative P-wave state with an inferred branching fraction [H] of B(B + — > ip(2s)K + — > J/iftr]K + ) = 
(2.1 ± 0.2) x 10" 5 . If we observe a two body S-wave 3/iftrj resonance, it would have the quantum 



S 



numbers, J PC = 1~ 1 , which are those of the unconfirmed h c . However, since the mass threshold 
of J/ipr] is 3.64 GeV/c 2 , which is higher than the lowest lying h c mass (estimated as the average 
of the three x state masses), such an observed resonance could indicate a new particle state. We 
note in addition that if a J/iprj resonance is found with a mass higher than open charm threshold 
(~3.77 GeV/c 2 ) it is unlikely to be a conventional charmonium state. 

The data used in this analysis correspond to a total integrated luminosity of 81.87 fb _1 taken 
on the T(45) resonance producing a sample of 89.96 ± 0.99 million BB events (N B jj). The data 
were collected at the PEP-II asymmetric-energy e + e~ storage ring with the BABAR detector, fully 
described elsewhere [7]. The BABAR detector includes a five-layer silicon vertex tracker (SVT) 
and a forty-layer drift chamber (DCH) in a 1.5-T solenoidal magnetic field. These devices detect 
charged particles and measure their momentum and energy loss. Photons and neutral hadrons 
are detected in a CsI(Tl) crystal electromagnetic calorimeter (EMC). The EMC detects photons 
with energies as low as 20 MeV and identifies electrons by their energy deposition. An internally 
reflecting ring-imaging Cherenkov detector (DIRC), composed of quartz bars, measures the charged 
particle velocity for particle identification. Penetrating muons and neutral hadrons are identified 
by the steel flux return, which is instrumented with 18-19 layers of resistive plate chambers (IFR). 

The preliminary selection criteria in this analysis follow previous BABAR analyses jH] and are 
briefly recalled here. All charged track candidates are required to have at least 12 DCH hits and 
transverse momentum greater than 100 MeV/c. The track candidates not associated with a K s 
decay must also originate near the nominal beam spot. The muon, electron, and kaon candidates 
must have a polar angle in radians of 0.3 < 6> M < 2.7, 0.410 < 6 e < 2.409, and 0.45 < 6 K < 2.50, 
respectively. In addition, all charged kaon candidates are required to have a laboratory momentum 
greater than 250 MeV/c. These requirements ensure the selection of tracks in the regions where 
the acceptance is well understood by the particle identification (PID) systems. 

The detailed explanation the PID is given elsewhere [5], jlOj . Briefly, photon candidates are 
identified from energy deposited in contiguous EMC crystals, summed together to form a cluster 
with total energy greater than 30 MeV and a shower shape consistent with that expected for 
electromagnetic showers. Electron candidates are required to have a good match between the 
expected and measured energy loss (dE/dx) in the SVT and DCH, and between the expected and 
measured Cherenkov angle in the DIRC. The measurements of the ratio of EMC shower energy to 
track momentum, and the number of EMC crystals associated with the track candidate must be 
appropriate for an electron. Muons are selected based on the energy deposited in the EMC, the 
number and distribution of hits in the IFR, the match between the IFR hits and the extrapolation 
of the track by the DCH into the IFR, and the depth of penetration of the track into the IFR. 
Charged kaon candidates are selected based on energy loss information from the SVT and DCH 
and the Cherenkov angle measured by the DIRC. 

The intermediate states in the decay modes used in this analysis, J/tp (e + e~), J/ip (/i + /i~), and 
Kg (7r + 7r~), are selected within the mass intervals listed in Tabled The di-electron mass interval 
is larger than the di-muon to account for Bremstrahlung in the detector. The K s decay length is 
required to be greater than 0.1 cm. 

The remaining four selection criteria are the rj mass, the ir° veto, the photon helicity angle [Hj 
from the r/ decay and the thrust angle. In the following we describe each of these criteria and the 
process used to determine their final cut values. 

The r\ candidates are required to have 77 mass within the range listed in Tabled If either of 
the photons associated with an r/ candidate, in combination with any other photon in the event, 
forms a 77 mass near the nominal 7r° mass, the rj candidate is vetoed as a tt° background. The 
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possible 77 masses and 
is near 1, where 



cos 0j£ 



7T° veto selection was adjusted by limiting the mass difference between a! 
the nominal ir° mass as shown Table ^ The r\ candidate is rejected if 
6^ is the photon helicity angle in the r\ rest frame. This rejects combinatoric background due to 
random pairs of photons that typically have a photon helicity angle distribution that peaks at or 
180 degrees. 

An additional requirement is applied to separate two-jet continuum events from the more spher- 
ical B meson decays. The angle 0t between the thrust |H] direction of the B meson candidate and 
the thrust direction of the remaining tracks in the event is calculated. We reject events when |cos 0t\ 
is near 1, since the distribution in cos^y is flat for BB events, while for background continuum 
events the distribution is peaked at cos 0t = ±1. 

Estimation of the signal and the background use two kinematic variables: the energy difference 
AE between the energy of the B candidate and the beam energy E% in the T(4S) rest frame; and 

the energy-substituted mass toes = "J (E£) 2 — (P B ) 2 , where P B is the reconstructed momentum 
of the B candidate in the T(45) frame. Typically these two variables in a two-dimensional plot 
for the B meson signal will appear as a weakly correlated two-dimensional Gaussian distribution, 
whereas the background is roughly uniformly distributed. The AE and tyies resolutions are mode 
dependent. A signal region for each mode, B + and B°, is defined as a rectangular region in the 
AE versus toes plane with |mES — msl < 7.5 MeV/c 2 , where tub is the mass of B meson and 
\AE\ <40 MeV. Before the data were analyzed, these four final selection criteria were optimized 
using a Monte Carlo (MC) simulation of the signal and the known backgrounds. Motivated by the 
B — > J /i/j(f)K measurement, the ab initio value of the branching fraction for B — ► J /ipr]K used 
in the signal MC was 5 x 10~ 5 . The number of reconstructed MC signal events (n™) and the 
number of reconstructed MC background events (n™ c ) in the signal box were obtained using the 
same criteria to estimate the signal significance ratio, n™ c / ^Jn™ c + n™ c . This ratio was maximized 
by varying these four selection criteria. The resulting criteria were fixed and are listed in Table [21 
The resulting AE and toes distributions for data are shown in Figs. Inland El The number of data 
events (no) observed in the signal box region for each mode, B + and B°, is listed in Table El 

The toes distribution background shape is determined by fitting the line shape of an ARGUS 
function |Hj for each mode to the mES distributions formed from the AE sideband region of the 
on-peak data of 0.1 < \AE\ < 0.14 for the B+ mode and 0.08 < \AE\ < 0.28 for the B° mode. To 
determine the number of signal events, we fit the data toes distribution to a model ( see Figs.[2fc) 
and EK C ) ) that is the sum of this ARGUS function whose normalization is allowed to vary, and 
a signal Gaussian, whose width is fixed to the MC value but whose mean and normalization are 
allowed to vary. The integral of this resulting ARGUS function over the signal box region yields 
the number of background events (n&) and its uncertainty (07,). Similar toes background shapes 
were found by using a 77 mass sideband outside the nominal r/ mass. Occasionally backgrounds 
with an inclusive J/ip decay will lead to a peaking background near or at the B mass in the ities 
distribution. Background shapes consisting of a sum of the ARGUS function and a broad Gaussian 
were attempted, however its contribution was found to be small and was dropped from the final 
background fit. The final values of rib and ab are listed in Table El 

We find evidence for signals in the J/tpr]K + and J/iprjKg modes as seen in Figs. El and El To 
determine the branching fraction for these modes the number of signal events (n s ) is given by a 
simple subtraction of the estimated number of background events from the events in the signal 
box region, n s = tiq — rib- The efficiencies (e) for each mode, listed in Table El are determined by 
MC simulation with three-body phase space and unpolarized J ftp decays. The calculation of the 
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Figure 2: The AE 1 and toes distributions for B + —* J/tpr]K + . The AE 1 vs. toes event distribution 
is shown in (a) with a small rectangle corresponding to the signal region selection defined in the 
text. The AE projection within a toes signal region selection is shown in (b). The toes projection 
within a AE signal region selection is shown in (c). The solid line in (c) is the fit described in the 
text. 
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Figure 3: The AE and toes distributions for B° — > J '/iprjKg. The descriptions of Figs. 3(a), (b) 
and (c) follow those of Figs. 2(a), (b) and (c), respectively. 



branching fraction {BF) is 

^ ^ IT'S 

BF = 

N BB X e X / 
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where e is the efficiency for the mode and / is the product of secondary branching fractions for the 
J/tp, 77, and Kg. The results are given in the last column of Table where the first and second BF 
uncertainities are the statistical and systematic, respectively. The statistical uncertainty is derived 
from the error on n s which is \JriQ + . 

In Table |1J we list the contribution to the systematic uncertainty from the error on each of 
the following quantities: Ng-^; secondary branching fractions jH]; MC statistics; PID, tracking, 
and photon detection efficiencies; 7r° veto; 77 mass range; background parameterization; and model 
dependence. The PID, tracking, and photon detection efficiency uncertainties are based on the 
study of data control samples The error in the 7r° veto efficiency was studied by measuring 
the inclusive 77 rate in both data and MC. After varying the 7r° veto mass cut, the change in the 
number of inclusive 77's was determined in data and MC to measure the deviation between the data 
and the simulation and to estimate the resulting systematic uncertainty. The uncertainty in the 
77 mass range was determined by comparing the measured 77 mass resolution in inclusive 77 decays 
to the 77 mass resolution from the signal MC. The background parameterization uncertainty was 
estimated by changing the ARGUS shape parameter by ±1 standard deviation, refitting the 7tt.es 
data distribution and recalculating the number of signal events. Additional systematic uncertain- 
ties due to the decay model dependence are estimated for the modes J/iprjK + and J/tjjrjKg. MC 
simulations are used to determine how much the efficiency depends on assumptions about interme- 
diate resonances and angular distributions. Five samples were studied. One sample is generated 
with 100% transversely polarized J/tp and another with 100% longitudinally polarized J/tp. The 
other three samples had large J/ipr] mass, large r]K mass or small J/ipK mass. The resulting 
relative change in efficiency is entered as a fractional systematic uncertainty in Tabled The total 
systematic uncertainty for each mode combines all these separate errors in quadrature, is listed the 
last column in Table 01 and is used to determine the BF systematic uncertainties listed in Table 01 

The probability for a null hypothesis (P-value) is defined as the Poisson probability that the 
estimated number of background events fluctuates to observed number of events uq or greater. In 
Table |3] we provide P-values calculated using the central value of the background estimate rib an d 
the value increased by 1 standard deviation, rib + crb, to provide an estimate of probability including 
the background systematic uncertainty. 

We determine also the 90% confidence level upper limit on the branching fraction using no, 
rib, an d Cft, in the signal region, and the total systematic uncertainty ax- Assuming the two 
uncertainties (07,, <tt) are uncorrelated and Gaussian, the Bayesian upper limit on the number 
of events (-/Vgo%) is obtained by folding the Poisson distribution with two normal distributions 
describing these two uncertainties and integrating the resulting function to the 90% confidence 
level (C.L.). This assumes that the a priori branching fraction distributions are uniform. The 
results are listed in Table 01 

Our preliminary branching fraction is comparable to B —* J /ip(j)K and is consistent with simple 
expectations from color-suppressed decay diagrams. The ratio of the charged (J/ijjrjK^ 1 ) to neutral 
(J/iprjKg) branching fractions is expected to be two which is not inconsistent with our results. 
Although the current statistics are too small, a larger data set of the neutral mode could be useful 
for a test of CP violation and both charged and neutral modes may enable a search for resonances 
in the two body intermediate states. 

In summary, we find evidence for the decay B — > J/ipr/K in two modes with preliminary 
branching fractions of, B(B + J/^r]K + ) = (10.8 ± 2.3 ± 2.4) x 10" 5 and B(B° -> J/ifyqK%) = 
(8.4 ±2.6 ±2.7) x 10" 5 . 
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TABLES 



Table 1: Mass regions for selection of intermediate particles. 



Mode Mass Range (GeV/c 2 ) 



J ftp — > e + e 


2.95 


< 


M(e+c-) 


< 


3.14 




3.06 


< 


M(m + M") 


< 


3.14 


K° s -► vr+vr- 


0.489 


< 


M(7T + 7T _ ) 


< 


0.507 


»7-+ 77 


0.525 


< 


M( 77 ) 


< 


0.571 



Table 2: Final selection criteria for the B 1 * 1 — > J /■ipr]('j^)K ± and B° — > J/V>r/(77)Kg modes. 



Variable 


J/4>T]K + 




\M V - 0.547| < 


0.023 GeV/c 2 


0.023 GeV/c 2 


7T° veto if |M( 7r , + lother ) - 0.135| < 


0.017 GeV/c 2 


0.010 GeV/c 2 


Helicity: \cos(9^)\ < 


0.93 


0.81 


Thrust: | cos (0 T )| < 


0.8 


0.9 



Table 3: Branching fractions and 90% C.L. upper limits. 



Mode 


e 


no 


n b ± o b 


^90% 


90% C.L.U.L. 

(io- s ) 




P-value 
range 


Branching Fraction 

(io- s ) 




10.75% 
8.53% 


99 
39 


50.3 ±3.0 
18.5 ± 1.7 


70.0 
34.5 


<15.5 
<14.1 


(0.09 
(.23 


- 1.42) x 10~ 8 

- 1.3) x 10~ 4 


10.81 ±2.31 ±2.37 
8.35 ±2.64 ±2.67 



Table 4: Systematic uncertainty summary on the branching fractions. All are fractional uncertain- 
ties in percent. 



Mode 


^BB 


Secondary 


M.C. 


PID, 


7T 


V 


Bkgd. 


Model 


Total 




Branching 


Stat. 


tracking, 


veto 


mass 


Param. 










Fraction 




7 Det. 




cut 










1.1 


2.48 


1.77 


8.2 


8.1 


3.40 


16.7 


5.1 


22.0 


J/i>r,K° s 


1.1 


2.52 


2.17 


8.3 


8.3 


3.14 


27.0 


9.5 


32.0 
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